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1 Introduction 



A well-known class of gravitational waves are the plane- fronted waves with par- 
allel rays ( pp- waves) which admit a covariantly constant null vector M , H , [0] , Q . 
The metric can be written using cylindrical coordinates: 

ds 2 = 2dudW + 2Hdu 2 - dP 2 - P 2 d<$> 2 (1) 

where H = H(u,P,&) and is of Petrov type N or conformally flat. Mainly 
axisymmetric pp-waves have been discussed in the literature. Such is the 
Schwarzschild solution boosted to the speed of light, which is interpreted as 
a massless null particle || , || or an ultrarelativistic black hole @ , || , |§ . A ring 



of massless particles is produced by boosting the Kerr metric [ 10 [ . Axisym- 



metric pp- waves describe also the gravitational field of light beams |ll|,[12|. 



Plane-fronted electromagnetic waves generate pp-gravitational waves as exact 



solutions of the Einstein-Maxwell equations [13], [14], [15]. 

While superposition of pp-waves running in the same direction is trivial, 
their collisions have been studied mainly for the subclass of plane waves [16 [,[17] 



in which the physical invariants are constant over the wave surfaces. The reason 
is that (1) is unsuitable to describe two approaching waves, but in the case of 



plane waves a transformation due to Rosen [18] converts (1) into the Szekeres 



line element |],(T|],[||] which can encompass two approaching waves and the 
region of their interaction. For a wave of constant polarization it reads 

ds 2 = Idudv - e~ u (e v dr 2 + e~ v dy 2 ) (2) 

where U and V are functions of the null coordinate u . Having in mind that 
\/2u = t — z, y/2v = t + z it is clear that (2) is a diagonal metric. Unfortunately, 
such waves have infinite extent and energy. The waves given by (1) with a 
suitably chosen H are finite in extent and energy but are not asymptotically 
flat and contain metric discontinuities for impulsive and shock waves. 

In the present paper we generalize the Rosen transformation to axisymmet- 
ric pp- waves and find for them a diagonal and asymptotically flat form like (2). 
This may be considered both as an alternative description and as a preparatory 
step before the investigation of their head-on collisions. 



2 Generalized Rosen transformation 

Let us investigate upon what conditions on H(u, P, the metric (1) can be 
diagonalized. Following Rosen, we do not change u at all. The requirement 
g uv = 1 , like in (2), is ensured by W = v + W\ (u,r,(p) . Next, the vanishing 
of g vv takes place when P v — ^> v — . The non-diagonal term g r( p disappears if 

P r P^ + P 2 $ r $^ = (3) 
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Obviously P r / , / and the simplest solution of (3) is P^ = <3? r = . 
Then the remaining non-diagonal terms disappear when 

W lr = P u P r W lip = P 2 <5> u % (4) 

with P (u,r) and (u, if) . Integrating the first equation in (4) and plugging it 
into the second we find a l.h.s. independent of r , unlike the r.h.s. The simplest 
solution is <& M = and = ip which does not change the range of the angular 
coordinate. Hence W\ = W\ (u, r) . The nullification of g uu is guaranteed by 
the relation 

2H = -2W lu + Pi (5) 

The r.h.s. depends on u and r , therefore after the coordinate transformation 
H = H(u,r), which means that in the beginning it was H(u,P) . Thus the 
sufficient condition for the diagonalization of (1) is the axial symmetry of H 
. We can't say that it is also necessary because the simplest solutions of (3,4) 
have been used. 

Integrating eq(4), inserting the result into (5) and taking the r-derivative we 
obtain the main equation which governs the generalized Rosen transformation: 

P(u,r) uu = -H(u,P) P (6) 

The line element (1) becomes 

ds 2 = 2dudv - P 2 dr 2 - P 2 d(p 2 (7) 

which is diagonal and the coordinate transformation reads 

u = u <$> = (p P = P(u,r) W = v + Wi(u,r) (8) 

where P is determined by (6) and W\ by (4). W does not appear in (6,7). 
Usually (1) is written in cartesian coordinates: 

ds 2 = 2dudW + 2H (u, X, Y) du 2 - dX 2 - dY 2 (9) 

For comparison, we apply the same diagonalization process to (9) with the 
following results. H must be separable 

H = H 1 (u,X)+H 2 (u,Y) (10) 

there are two main equations 

X (u, x) uu = -H, (n, X) x Y (n, y) uu = -H 2 (u, Y) Y (11) 

the line element becomes 

ds 2 = 2dudv - X 2 dx 2 - Y 2 dy 2 (12) 

and the coordinate transformation reads 

u = u X = X(u,x) Y = Y(u,y) 
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W = v + J X X u X x ,dx' + J V Y u Y y ,dy' (13) 

In vacuum the only surviving Einstein equation gives 

H 1X x + H 2YY = (14) 

After the usual removal of linear terms in Hi (11) becomes linear: 

X uu = -f(u)X Y = f{u)Y (15) 

where / (u) is the second derivative of H% . This permits the choice X = xF (u), 
Y = yG (u) which represents exactly the Rosen transformation for plane waves 
of constant polarization pjl, [pq1 ,pC|[. Thus in vacuum the diagonalization of 
pp-w&ves in cartesian coordinates requires separability (10) and leads directly 
to plane waves. 



3 pp-waves in Brinkmann and diagonal form 

We have shown that an axisymmetric pp-w&ve can be described by (7). At first 
sight (7) is a special case of the Szekeres line element (2) where now U and V 
depend on u and r. In fact, there is no loss of generality and any metric (2) 
may be written as (7) provided the Einstein equations hold. Let us compare 
these equations for (1,7) and (2). 

A pp-wave allows energy-momentum tensors of few types: vacuum, null 
electromagnetic field or pure radiation (null dust), which may be combined 
01, |H . All of them have only one non-trivial component: 

T uu = 2p{u,P) = 2p R + 2 PE (16) 

where pe is the energy-density of pure radiation with no matter equations and 
Pe is the electromagnetic energy-density 

2 PE = V^VV if> = A u (17) 

We have used relativistic units with 8-7rG/c 4 = 1 . A u is the only component of 
the vector potential and satisfies the Maxwell equation Aip = . We suppose 
that no charges and currents are present. The gradient and Laplacian are with 
respect to P, $ (or X, Y ). We have chosen this formalism instead of the 
Newman-Penrose one since T uu , A u and the Maxwell equation do not change 
under the generalized Rosen transformation (8). 

In Brinkmann coordinates the only non-trivial Ricci tensor component is 
R uu which is a Laplacian. The Einstein equation is 

Hpp + ^H P = 2p (18) 
and its solution is well-known from classical potential theory: 

R P = \\l P («' P ') P ' dP ' + J>P* i U ) ( 19 ) 
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f p d p> r p> p 

H = 2 -5T P ") P " dP " + 2 P- ( u ) ln - ( 2 °) 

where an ignorable term has been omitted in (20). The second term in (19,20) 
is the exterior solution with arbitrary p e and some constant length a . The 
Maxwell equation is 

(PiPp)p + pVw = ( 21 ) 

We have retained some ^-dependence in tj} but it must disappear in pe ■ The 
only non-zero Weyl scalar is given by 

Hp 

*4 = P-^ (22) 

Now, let us concentrate on the metric (2) with U (u, r) , V (u, r) . The 
Einstein equations yield: 

2U UU = Ul + Vl + 4p (23) 

(U + V) ru = (U + V) r V u (24) 
(U + V) rr = (U + V) r V r (25) 
Eqs(24,25) are easily integrated 

(U + V) r = -2e v (26) 

We have chosen the integration constant in order to restore Minkowski space- 
time for V = U , r = e~ v . This allows smooth transition to it in front of the 
wave and is in accord with our demand for asymptotically flat solutions. The 
Maxwell equation is 

(e- u i/Jr) + e^VW = ( 27 ) 

where like in (21) we allow for some (/^-dependence. Using the natural NP null 
tetrad f§ we find three non-trivial Weyl scalars: 

*2 = -& u - v [(U + V) rr -(U + V) r V r ] * 3 = -^e^Rur (28) 

^4 = \ (VuUu - Vuu) (29) 

A look at the Einstein equations (23-25) shows that ^ = ^3 = and the field 
is of type N not of type II. 

One can try to solve the relevant eqs(23,26), when p is given , in two ways. 
First, we may take an arbitrary U , solve for V from (23) and insert the result 
into (26). This leads to a condition on U : 

(U + A) r = ai (r) e A + a 2 (r) A = J ^2U UU - U* - Apdu (30) 

with arbitrary a± / and 02 . Eq(30) is too complicated to be examined. 
Second, we take an arbitrary V and notice that (23) is a linear second-order 
equation for e~ u ^ 2 . Then we integrate (26): 

U = - ( (2e v + V r ) dr + f x (u) (31) 
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where f\ (u) is some yet undetermined function. Substituting (31) into (23) we 
get an equation for /i with additional conditions on V to yield a r-independent 
fi which again are very complicated. 

Let us now compare the two theories. The link is given by (7): 

P = e~^ P r = e^ (32) 

Eq(32) gives at first sight an additional constraint between U and V but this 
turns out to be exactly (26) which is necessarily satisfied. Going backwards, 
(26) shows that (32) holds for some P . It can be shown further that (6,18) 
are equivalent via (32) to (23). The same is true about (21) and (27) which 
is not so surprising for a Laplacian. At last, under (32) the Weyl scalar (22) 
coincides with (29). Consequently, axisymmetric pp-waves (1) are in a one-to- 
one correspondence with the solutions U (u, r), V (u, r) for metric (2). Thus we 
can replace metric (2) with two functions by metric (7) with one function P or 
by metric (1) with one function H . Each of these forms has its own merits. 
The Brinkmann metric (1) has simple Einstein equations (compare (18) with 
(30,31)) but is not asymptotically flat for exterior solutions, sometimes has a 
discontinuous H and is unfit for studying collisions of pp- waves because of only 
one null coordinate. The metric given by (7) or (2,32) is worth as a starting 
point for the interaction problem and is asymptotically flat for realistic p as will 
be shown in the following. However, < P < oo as a radial coordinate in (1) 
which makes g w in (7) singular at some points. This coordinate singularity is 
innocuous when it is due to the cylindrical character of the coordinate system. 
If not, the experience with plane waves teaches that it becomes a fold singularity 
and is intimately related to the curvature singularity in the interaction region 



[21]. 



4 Solutions: general features 

Eq(6) with H satisfying (20) is a second-order nonlinear differential equation 
with respect to P . In the process of solving it arbitrary functions of r arise 
which reflect the residual freedom in the coordinate transformation and may be 
selected to further simplify the solution and satisfy boundary conditions. 

The trivial Minkowski solution is given hyH = 0,P = r,P r = l. Having 
in mind the setting of the collision problem, u = must be the boundary 
between the running wave ( u > ) and Minkowski spacetime ( u < ) where 
the wave has not yet arrived. This gives the universal boundary condition 

P(0,r)=r (33) 

We also demand asymptotic flatness i.e. P (u, oo) = r for fixed u . 

It is clear from (20) that the exterior solution is always separable, the interior 
is separable when p = p\ (P) p2 (u) ■ Almost all pp- waves discussed in the 
literature are of this kind with H (u, P) = H\ (P) p2 (u) and we shall consider 
only them in the following. A natural question arises: when H is static in (1) 
is it possible that P is also static in (7)? This is not allowed by (4,5). A static 
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P has P u = , W\ = W\ (u) and H = H (it) contrary to our assumption that 
H = H (P (r)) . Hence P depends on u even when p 2 (u) = 1 . 
For the exterior solution given in (19) eqs(6,22) become 

Puu = ~j>Pe («) ^4 = -jnPe («) (34) 

Unlike (15), eq(34) is non-linear and we can't get rid of the r-dependence. For 
many simple choices of p e (34) falls in the class of Emden- Fowler equations |j22|| . 
They are quite difficult to solve and many of them remain non-integrable. For 
example, when p e (u) = u n the integrable cases are just n = 0;— 1;— 2 . For 
asymptotically flat solutions (34) shows that ^4 — > when r — ► 00 . 
The general separable interior solution emerges from 

P U u = -H 1 (P) P p 2 (u) (35) 

and again is reducible in many cases to the Emden-Fowler equations and their 
generalizations. The case p\ (P) = 1 is special. Then (35) is linear in P and we 
can use an analog of the ansatz applied after (15), namely P (u,r) = rp(u) : 

H = \ (X 2 + Y 2 ) p 2 (u) ^4 = (36) 

Puu = -P2(u)p (37) 

ds 2 = 2dudv - p {uf (dX 2 + eZY 2 ) (38) 

This is the case of pure radiation with density pu = p 2 or an electromagnetic 
wave with potential, Maxwell scalar and energy-density given by 

if) = a 3 (it) X + a 4 («) Y $ 2 («) = - [03 («) - ia 4 («)] (39) 

PE (u) = p 2 («) = - (o§ + al) (40) 

where 03 , 04 are arbitrary functions. Obviously ip depends on (p while pe does 
not. The potential satisfies the Maxwell equation (21). In fact eqs(36-40) rep- 
resent a plane electromagnetic wave [|| and an axisymmetric electromagnetic 
pp-w&ve at the same time. There is no pure gravitational wave in addition 
because ^4 = 0. The Ricci scalar $22 = $2^2 is constant over the wave sur- 
face. On the contrary, pure plane gravitational waves can not be axisymmetric 
because their H ~ X 2 — Y 2 which is ^-dependent. 

The case discussed above provides a link between plane and axisymmetric 
pp-waves. Even for it, eq(37) is the normal form of the general linear second- 
order equation and its general solution is given analytically only if a non-trivial 
concrete solution is known. Therefore we are going to discuss two cases of 
simple li-dependence when the solution of (6) may be found. These are the 
impulsive and shock waves. 
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5 Impulsive waves 



These are waves with p2 = 5 (u) . Eq(37) may be integrated with the help of 
(33): 

P = r(l- Stu) P r = 1 - H lrr u (41) 

Eq(22) transforms into 



l Fi 



Pi ( r ) ~ — ) 5 (u) 



(42) 



which clearly demonstrates the impulsive character of the wave. It is seen 
from (19) that H\ r > and (41) shows that P always possesses a coordinate 
singularity for some u > 0, different from the cylindrical singularity at r = 
. This is a consequence of the positive energy condition and the idealized 
impulsive character of the wave. 

For a boosted Schwarzschild solution || , || , || H — 2/j.S (u) InP 2 where \i is 
the momentum of the null point-like particle and (41,42) give 



ds = 2dudv 



1 + ^u6 («) 



dr z 



1 



^u9 (u) 



r 2 dy 2 



* i 



8(r) 



fj,S (u) 



(43) 



(44) 



Eq(43) is exactly the line element found in @,||. There is a curvature singu- 
larity at the point of the source r = . H is also the function for an exterior 
impulsive solution given in (20). If t is fixed, for any z and r — > oo the solution 
is asymptotically flat. There is a coordinate singularity at \^2r 2 = 4/x(t — z) . 
For a fixed z , as time goes by, the singular circle centred at z expands towards 
infinity. 

For a boosted Kerr solution [lCtl : 



H = 2fi5 («) In 



According to (41): 

P : 



(i - 



Pr 



(45) 



(46) 



{r 2 -b 2 Y 

where b is the radius of the ring of massless particles. The curvature singularity 
moves to r = b and the region r < b is free of coordinate singularities. The 
metric is asymptotically flat. 

As a final example we present the diagonalization of an impulsive beam of 



light with transverse radius a |ll|, p2j| . This is a global solution the interior 
being given by (36) and the exterior by (20). The junction conditions require 
that 

Am P 2 f P \ 

H = —9 (a - P) 5 (u) + Am l + 21n- ) (P - a) 5 (u) (47) 



a* \ a J 

where m is the constant energy density. With the help of (18) eq(41) may be 
rewritten as 



(l-^t 



Pr 



1 + ( % 



2/0i u 



(48) 
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Inserting (47) into (48) we obtain for the interior and exterior solutions: 



P 



(l " » 
(l - ^u) 



p. 



8m 



U 



1 + 



8m 



U 



(49) 
(50) 



It is seen that P is continuous at r = a but P r makes a finite jump. Accord- 
ing to (48) the reason is the jump in p\ from zero to a finite constant, since 
the junction conditions require that Hi and H\ r should be continuous. Con- 
sequently, solutions which are perfectly well joined in Brinkmann coordinates 
acquire discontinuous metric upon diagonalization due to unrealistic densities 
with 6 (r) terms. The problem disappears when the density smoothly falls to 
zero. Take for example p\ (P) = e~ p . Then 



H 



InP 



S(u) 



(51) 



1 



Pr 



1 + 



-2r 2 e~ 



-U 



(52) 



When P — ► 0,oo H in (51) approaches the first or the second term in (47). 
Correspondingly, when r — > (52) approaches (49) and when r — > oo it ap- 
proaches (50) with 8m = a = 1 . The metric (52) is asymptotically flat but the 
coordinate singularities still exist. 



6 Shock waves 

These waves have H = H\ (P) 9 (u) and (6) has a first integral: 

Pl = c(r)-2H 1 (P) (53) 

It is clear from (20) that Hi is a positive and increasing function. This is the 
reason to keep the arbitrary function c (r) in (53) so that the r.h.s. is positive. 
Eq(53) is easily integrated. Imposing (33) we obtain 

f p dP' 

±u =i. 7m=mm =K{p - r) - K(r) (54) 

which gives P (u,r) indirectly. For future convenience we have introduced also 
the indefinite integral K . 

In order to understand the meaning of c (r) let us discuss the interior solution 
(36-38) with p2 (u) = 6 (u) = 1 in the region occupied by the wave. The integral 
in (54) can be evaluated: 

P r 
arcsin — — = —u + arcsin —= (55) 
VC Vc 

and P is found by inverting (55). Let us choose 

c (r) = 2Hi (r) (56) 
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Then we obtain 



r cos u 



ds z 



2dudv — cos 2 u (dr 2 + r 2 dip z 



(57) 
(58) 



This, however, is the line element of an electromagnetic shock wave with Ricci 
scalar $22 = (u) Q and this is really the case here because we can choose 
111 (39) a 3 = y/26 (u), 04 = 0, $2 = —@(u) ■ We conclude that plane waves 
recommend the receipt (56). Eq(57) shows that it is equivalent to the method 
used in (35-37) for the linear case. 

Let us apply this receipt to the exterior solution H\ (P) = b In , b > . 
The integral in (54) yields the error function 



erf \ / In 




This formula may be inverted 



r exp 



P 



erf" 1 W 



(59) 




(60) 



The metric satisfies the necessary boundary condition (33) and is asymptotically 
flat for fixed u . The problem is that erfz < 1 and (59) imposes the constraint 



u < 




(61) 



The solution (60) does not cover the whole region u>0,0<r<oo. 

This is a generic feature of the choice (56) into (54). Then P < r because H\ 
is an increasing function. Hence, the minus sign must be chosen in (54). When 
u increases P necessarily decreases, becomes null and sometimes even negative 
as (57) demonstrates. However, for fixed r it remains bounded in order to keep 
the root in (54) real. The integral in (54) also remains finite so there should 
be some limit for the growth of u like (61). The same happens in (55) if we 
stick to the main branch of arcsin x . Fortunately, u (P) may be a multivalued 
function while P (u) can not be. This explains why there are no problems in this 
case. Multivalued functions resulted from the inversion of periodic functions. 
In the general case periodic functions do not appear in K and that causes the 
limit problem. If P is not extended to negative values the limit of u is also a 
coordinate singularity and is given by 



u 



K(r)-K(0) 



(62) 



This singularity is present generically in solutions with (56). 
Another choice is 

c(r) = 2H = 2H 1 (P ) (63) 

where Hq is some very big constant. Then we may take the positive sign in 
(54) and Pq > P > r . The P , u and r dependencies separate: 



K(P) = u + K(r) 



(64) 
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There is no coordinate singularity but the region P > Pq is not described by 
this coordinate system because K is ill-defined there. In turn this means 

r <P u<K(P )-K (0) (65) 

For the exterior solution these inequalities look like 

r < ae~b~ u < ay^e^erf <J ^ — In £ (66) 

If the first of them is made stronger and Hq is taken big enough, u and r can 
cover a lot of their range. The choice (63) is perfect if H\ (P) were bounded 
from above and Hq > ffi max • Unfortunately, this does not happen due to 
the lower limit of the inside integral in (20). It cures the singular behaviour at 
small P but generates a logarithmic term in H\ like the first term in (51). 

A problem arises when we try to join the interior and exterior solutions 
discussed above. We start with (47) and 5 (u) replaced by 9 (u) pT|| . Now we 
can't replace e.g. 9 (a — P) by 9 (a — r) and that makes eq(54) intractable. Like 
in the case of impulsive waves it is preferable to have one smoothly falling out 
p for all r like the example given by (51). With such Hi the integral in (54) 
can not be evaluated analytically and the limit problem still exists. This is the 
best we can do for realistic shock waves. 

Acknowledgement 

This work was supported by the Bulgarian National Fund for Scientific Research 
under contract F-632. 



11 



References 

[1] Brinkmann M W 1923 Proc. Nat. Acad. Sci. U.S.A. 9 1 
[2] Peres A 1959 Phys. Rev. Lett. 3 571 

[3] Ehlers J and Kundt W 1962 in Gravitation: an Introduction to Current 
research (ed.L.Witten, New York: Wiley) 

[4] Griffiths J B 1991 Colliding Plane Waves in General Relativity (Oxford: 
Clarendon) 

[5] Aichelburg P C and Sexl R U 1971 Gen. Rel. Grav. 2 303 

[6] Dray T and 't Hooft G 1985 Nucl. Phys. B 253 173 

[7] D'Eath P D 1978 Phys. Rev. D 18 990 

[8] D'Eath P D and Payne P N 1993 Phys. Rev. D 46 658 

[9] Bonnor W B, Griffiths J B and MacCallum M A H 1994 Gen. Rel. Grav. 
26 687 

[10] Ferrari V and Pendenza P 1990 Gen. Rel. Grav. 22 1105 
[11] Bonnor W B 1969 Comm. Math. Phys. 13 163 

[12] Ferrari V, Pendenza P and Veneziano G 1988 Gen. Rel. Grav. 20 1185 

[13] Peres A 1960 Phys. Rev. 118 1105 

[14] Bonnor W B 1970 Int. J. Theor. Phys. 3 57 

[15] Kramer D, Stephani H, Herlt E and MacCallum M 1980 Exact Solutions 
of Einstein's Field Equations (Cambridge: Cambridge University Press) 

[16] Baldwin O R and Jeffery G B 1926 Proc. Roy. Soc. A 111 95 

[17] Bondi H and Pirani F A E 1959 Proc. Roy. Soc. A 251 519 

[18] Rosen N 1937 Phys. Z. Sowjet. 12 366 

[19] Szekeres P 1970 Nature 228 1183 

[20] Khan K A and Penrose R 1971 Nature 229 185 

[21] Matzner R A and Tipler F J 1984 Phys. Rev. D 29 1575 

[22] Polyanin A D and Zaitsev V F 1995 Handbook of Exact Solutions for Or- 
dinary Differential Equations (Boca Raton: CRC Press) 



12 



